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Abstract 

Cyclic liftings are proposed to lower the error floor of low-density parity-check (LDPC) codes. The liftings 
£f) • are designed to eliminate dominant trapping sets of the base code by removing the short cycles which form the 
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a cycle c of length £(c) in the base graph such that the inverse image of c in the lifted graph consists of only 
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cycles of length strictly larger than 1(c). The proposed method is universal in the sense that it can be applied to 
any LDPC code over any channel and for any iterative decoding algorithm. It also preserves important properties 
of the base code such as degree distributions, encoder and decoder structure, and in some cases, the code rate. The 
proposed method is applied to both structured and random codes over the binary symmetric channel (BSC). The 
error floor improves consistently by increasing the lifting degree, and the results show significant improvements in 
the error floor compared to the base code, a random code of the same degree distribution and block length, and a 
random lifting of the same degree. Similar improvements are also observed when the codes designed for the BSC 
are applied to the additive white Gaussian noise (AWGN) channel. 
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I. Introduction 

Low-density parity-check (LDPC) codes [H have emerged as one of the top contenders for capacity 
approaching error correction over many important channels. They not only perform superbly but also 
lend themselves well to highly efficient parallel decoding algorithms. A well-known construction of 
LDPC codes is based on protographs, also referred to as base graphs or projected graphs [15J. In such 
constructions, a bipartite base graph G is copied N times and for each edge e of G, a permutation is 
applied to the N copies of e to interconnect the N copies of G. The resulting graph, called the N -cover 
or the N -lifting of G, is then used as the Tanner graph IfTTTl of the LDPC code. If the permutations are 
cyclic, the resulting LDPC code is called quasi-cyclic (QC). QC LDPC codes are attractive due to their 
simple implementation and analysis lfT5ll . 

At very large block lengths, the performance of LDPC codes can be well estimated using asymptotic 
techniques such as density evolution [14J. At finite lengths, however, our understanding of the dynamics of 
iterative decoding algorithms is limited. In particular, iteratively decoded finite-length codes demonstrate 
an abrupt change in their error rate curves, referred to as error floor, in the high signal to noise ratio (SNR) 
region. The analysis of the error floor and techniques to improve the error floor performance of LDPC 
codes are still very active areas of research. For the binary erasure channel (BEC), the error floor is well 
understood and is known to be caused by graphical structures called stopping sets flU. Richardson related 
the error rate performance of LDPC codes on the binary symmetric channel (BSC) and the additive white 
Gaussian noise (AWGN) channel to more general graphical structures, called trapping sets, and devised 
a technique to estimate the error floor lfT3l . Other estimation techniques based on finding the dominant 
trapping sets were also proposed for the BSC in Q and for the AWGN channel in |BI. Ifl6l In E2l - 
Xiao and Banihashemi took a different approach, and instead of focusing on trapping sets which are the 
eventual result of the decoder failure, focussed on the input error patterns that cause the decoder to fail. 
A simple technique for estimating the frame error rate (FER) and the bit error rate (BER) of finite-length 
LDPC codes over the BSC was developed in Il22ll . The complexity of this algorithm was then reduced 
in Il24ll . and the estimation technique was extended to the AWGN channel with quantized output in [|23l . 
More recent work on the estimation of the error floor of LDPC codes is presented in [0, 0, to which 
the reader is also referred for a more comprehensive list of references. 

There is extensive literature on reducing the error floor of finite-length LDPC codes over different 
channels and for different iterative decoding algorithms. One category of such literature, focusses on 
modification of iterative decoding algorithms, see, e.g., Q, while another category is concerned with the 
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code construction. In the second category, some researchers use indirect measures such as girth [fT8l or 
approximate cycle extrinsic message degree (ACE) ||T9ll , while others work with direct measures of error 
floor performance such as the distribution of stopping sets or trapping sets li20l . IfTOl . [fTTfl. In [20], edge 
swapping is proposed as a technique to increase the stopping distance of an LDPC code, and thus to 
improve its error floor performance over the BEC. Random cyclic liftings are also studied in EDI and 
shown to improve the average performance of the ensemble in the error floor region compared to the base 
code. Ivkovic et al. IfTOl apply the same technique of edge swapping between two copies of a base LDPC 
code to eliminate the dominant trapping sets of the base code over the BSC. 

In the approach proposed here also, we focus on dominant trapping sets which are the main contributors 
to the error floor. We start from the code whose error floor is to be improved, as the base code. We then 
construct a new code by cyclically lifting the base code. The lifting is designed carefully to eliminate 
the dominant trapping sets of the base graph. This is achieved by removing the short cycles which form 
the dominant trapping sets. Our work has similarities to IfTOl and Il20l . The similarity with both IfTOl 
and 11201 is that we also use graph covers or liftings to improve the error floor performance of a base 
code. It however differs from IfTOl in that we restrict ourselves to cyclic liftings that are advantageous in 
implementation. Moreover, to eliminate the dominant trapping sets we use a different approach than the 
one in IfTOl . More specifically, our approach is based on the elimination of the short cycles involved in the 
trapping sets. To do so, we derive a necessary and sufficient condition for the problematic cycles of the 
base code such that they are mapped to strictly larger cycles in the lifted code. The difference with ll20l 
is that while ||20l is focused on the ensemble performance of random liftings, our work is concerned with 
the intentional design of a particular cyclic lifting. 

Given a base code and its dominant trapping sets over a certain channel and under a specific iterative 
decoding algorithm, the proposed construction can lower the error floor by increasing the block length 
while preserving the important properties of the base code such as degree distributions, and the encoder 
and decoder structure. The code rate is also preserved or is decreased slightly depending on the rank 
deficiency of the parity-check matrix of the base code. Moreover, the cyclic nature of the lifting makes it 
implementation friendly. We apply the proposed construction to the Tanner code |fT8l and two randomly 
constructed codes, one regular and the other irregular, to improve the error floor performance of Gallager 
A/B algorithms over the BScQ Simulation results show a consistent improvement in the error floor 

'The choice of BSC/Gallager algorithms is for simplicity, and the proposed construction is applicable to any channel/decoding algorithm 
combination as long as the dominant trapping sets are known. 
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performance by increasing the degree of liftings. The constructed codes are far superior to similar random 
codes or codes constructed by random liftings in the error floor region. We also examine the performance 
of the codes constructed for BSC/Gallager B algorithm, over the AWGN channel with min-sum decoding 
and observe similar improvements in the error floor performance. 

The remaining of the paper is organized as follows: Section UU introduces definitions, notations and 
background material used throughout the paper. In Section III, the proposed construction is explained and 
discussed. Numerical results are presented in Section [IV] and finally Section |V] concludes the paper. 

II. PRELIMINARIES: LDPC CODES, TANNER GRAPHS, GRAPH LIFTINGS AND 

TRAPPING SETS 

A. LDPC Codes and Tanner Graphs 

Consider a binary LDPC code C represented by a Tanner graph G = (VbUV c ,E), where Vb = 
{&i,..., b n } and V c = {ci, . . . , c m } are the sets of variable nodes and check nodes, respectively, and 
E is the set of edges. Corresponding to G, we have an m x n parity-check matrix H = [hij] of C, where 
hij = 1 if and only if (iff) the node Cj E V c is connected to the node bj E Vb in G; or equivalently, iff 
{bj, Ci} E E. If all the nodes in the set V& have the same degree d v and all the nodes in the set V c have 
the same degree d c , the corresponding LDPC code is called a regular (d v , d c ) code. Otherwise, it is called 
irregular. 

A subgraph of G is a path of length k if it consists of a sequence of k + 1 nodes {u±, . . . , Uk+i} and 
k distinct edges {{wj, : % = 1, . . . , k}. We say two nodes are connected if there is a path between 
them. A path is a cycle if u\ = Uk+i, and all the other nodes are distinct. The length of the shortest 
cycle(s) in the graph is called girth. In bipartite graphs, including Tanner graphs, all cycles have even 
lengths. So, the girth is an even number. 

B. Graph Liftings 

Consider the set of all possible permutations Sn over the set of integer numbers — {1, . . . , N}. 

This set forms a group, known as the symmetric group, under composition. Each element n E Sn can 
be represented by all the values i E Z\^n- For the identity element ir , we have 7To(i) = i,Vz, and 
the inverse of n is denoted by 7r _1 and defined as 7r _1 (7r) = it . It is easy to see that the symmetric 
group is not Abelian. An alternate representation of permutations is to represent a permutation ix with an 
N x N matrix II = [7^], whose elements are defined by Wij — 1 if j — 7r(i), and -n^ = 0, otherwise. 
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As a result, we have the isomorphic group of all N x N permutation matrices with the group operation 
defined as matrix multiplication, and the identity element equal to the identity matrix I N . In particular, 
corresponding to the composition n'(ir), we have the matrix multiplication II x II'. Moreover, since the 
permutation matrices are orthogonal, the inverse of a permutation matrix is its transpose, i.e., II -1 = II T . 

Consider the cyclic subgroup Cn of Sn consisting of the N circulant permutations defined by n d (i) = 
i + d mod N + 1, d G Z _».jv-i- The permutation n d corresponds to d cyclic shifts to the right. In the 
matrix representation, permutation ltd corresponds to a permutation matrix whose rows are obtained by 
cyclically shifting all the rows of the identity matrix I N by d to the right. This matrix is denoted by 1^ . 
Note that = I N . For the composition of permutations, we have 

j(di) x j(d 2 ) _ j{d\-\~d2 mod N) (j) 

Clearly, C N can be generated by where each element I^ d \ d G Z ^ N _i of C N is the d-th power of 
This defines a natural isomorphism between Cn and the set of integers modulo TV, Z Q ^ N -i, under 
addition. The latter group is denoted by Z N . 

Consider the following construction of a graph G = (V, E) from a graph G = (V, E): We first make 
N copies of G such that for each node v G V, we have N copies v = {v i . . . , v N } in V. For each edge 
e = {u,v} E E, we assign a permutation 7r e e S'at to the copies of e in E such that an edge {u i: Vj} 
belongs to E iff 7r e (i) = j. The set of these edges is denoted by e. The graph G is called an N -cover or 
an N -lifting of G, and G is referred to as the base graph, protograph or projected graph corresponding 
to G. We also call the application of a permutation 7r e to the N copies of e, edge swapping, high lighting 
the fact that the permutation swaps edges among the N copies of the base graph. 

In this work, G is a Tanner graph, and we define the edge permutations from the variable side to 
the check side, i.e., the set of edges e in E corresponding to an edge e = {b,c} e E are defined by 
{h, Ctt^)}, i G Zi->n- Equivalently, e can be described by {fe^-i^), Cj}, j G ^i^at. Our focus in this 
paper is on cyclic liftings of G, where the edge permutations are selected from C N , or equivalently Z N . 
Thus the nomenclature cyclic edge swapping. In this case, if 1^ is a permutation matrix from variable 
nodes to check nodes, I^; d! = N — d mod N, will be the corresponding permutation matrix from 
check nodes to variable nodes. It is important to distinguish between the two cases when we compose 
permutations on a directed path. 

To the lifted graph G, we associate an LDPC code C, referred to as the lifted code, such that the 
mN x nN parity-check matrix H of C is equal to the adjacency matrix of G. More specifically, H 
consists of m x n sub-matrices [H]ij , 1 < % < m, 1 < j < n, arranged in m rows and n columns. The 



5 



sub-matrix [H]ij in row % and column j is the permutation matrix from Cn corresponding to the edge 
{bj, Ci} where hij ^ 0; otherwise, [H]ij is the all-zero matrix. Let the m x n matrix D = [dy] be defined 
by [H]ij = I^', dij G Zn if ^ 0, and dij = +oo, otherwise. Matrix D, called the matrix of edge 
permutation indices, fully describes H and thus the cyclically lifted code C. 

C. Trapping Sets and Error Floor 

It is well-known that trapping sets are the culprits in the error floor region of iterative decoding 
algorithms. An (a, b) trapping set is defined as a set of a variable nodes which have b check nodes of 
odd degree in their induced subgraph. Among trapping sets, the most harmful ones are called dominant. 
Trapping sets depend not only on the Tanner graph of the code but also on the channel and the iterative 
decoding algorithm. In general, finding all the trapping sets is a hard problem, and one often needs to 
resort to efficient search techniques to obtain the dominant trapping sets [[T3l , 01, Il24ll . fl2T|. Trapping 
sets for Gallager A/B algorithms over the BSC are examined for a number of LDPC codes in 0, ll24l . 
In this work, we assume that the dominant trapping sets are known and available. 

In the context of symmetric decoders over the BSC, the error floor of FER can be estimated by 11221 

FER ^Nje J , (2) 

where e is the channel crossover probability, and J and Nj are the size and the number of the smallest 
error patterns that the decoder fails to correct. From ©, it is clear that the dominant trapping sets over 
the BSC are those caused by the minimum number of initial errors. (In 0, |[T0l . parameter J in © is 
called minimum critical number.) In the double-logarithmic plane, one can see from © that log(FER) 
decreases linearly with log(e) and the slope of the line is determined by J. 

III. DESIGN OF CYCLIC LIFTINGS TO ELIMINATE TRAPPING SETS 

In this work, our focus is on the design of cyclic liftings of a given Tanner graph to eliminate its 
dominant trapping sets with respect to a given channel/decoding algorithm with the purpose of reducing 
the error floor. (This is equivalent to the design of non-infinity edge permutation indices of matrix D.) 
For example, Equation © indicates that for improving the error floor on the BSC, one needs to increase 
J and decrease Nj corresponding to the dominant trapping sets. In particular, while increasing J for 
dominant trapping sets increases the slope of log(FER) vs. log(e) at low channel crossover probabilities 
e, reducing Nj amounts to a downward shift of the curve. So, the general idea is to primarily increase 
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J by eliminating the trapping sets with the smallest critical number, and the secondary goal is then to 
decrease Nj. 

It is well-known that dominant trapping sets composed of short cycles 0, Il24l . To eliminate the 
trapping sets, we thus aim at eliminating their constituent cycles in the lifted graph. In the following, we 
examine the inverse image of a (base) cycle in the cyclically lifted graph. 

A. Cyclic Liftings of Cycles 

Lemma 1: Consider a cyclic iV-lifting G of a Tanner graph G. Consider a path £ of length t in G, 
which starts from a variable node b and ends at a variable node b' with the sequence of edges e 1 , . . . ,e£. 
Corresponding to the edges, we have the sequence of permutation matrices I^ dl \ . . . , I^ dl >. Then the 
permutation matrix that maps b to b' through the path £ is I^ d \ where 



Proof: The permutation matrix that maps b to b' is obtained by multiplying the permutation matrices 



of the "directed" edges along the path. This results in 1^ = 1^ x i(N-d 2 modN) x . . . x j(iV-^modAr) 



The value of d given in ([3]) is called the permutation index of the path from b to b' . Clearly, the 
permutation index of the path from b' to b is equal to d! = N — d mod N. If b = b' and all the other 
nodes are distinct, then the path will become a cycle and depending on the direction of the cycle, its 
permutation index will be equal to d or d'. 

Theorem 1: Consider the cyclic iV-lifting G of the Tanner graph G. Suppose that c is a cycle of length 
£ in G. The inverse image of c in G is then the union of N/k cycles, each of length k£, where k is the 
order of the element(s) of Zn corresponding to the permutation indices of c. 

Proof: We first note that the elements of Z N corresponding to both permutation indices of c have 
the same order. Suppose that the permutation indices of c in the two directions are equal to d and 
d! = N — d mod N. Now, d! is the inverse of d in Z N , and thus has the same order as d. 

Denote the sequence of variable and check nodes in c by bi i: c i2 , 6j 3 , c i4 , . . . , Ci i( b ix . Starting from 
any of the iV variable nodes in b^, say 6j U , 1 < j < N, without loss of generality, we follow one of the 
two paths of length i in the inverse image of c in G, which corresponds to the direction on c associated 
with the permutation index d. As a result, based on Lemma [Q we will end up at the variable node b ilP 
in 6 il9 where p = j + d mod iV and d is given in ©. If d = 0, then the path ends at b ix j, meaning 




(3) 



The lemma is then proved using (OQ). 
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that the inverse image of c which passes through b ix j is a cycle of length £. Similarly one can see that 
the inverse image of c passing through any of the N nodes fe^i, . . . , b^N is a cycle of length £, and since 
these cycles do not overlap, the inverse image of c in this case is the union of N cycles, each of length 
£. This corresponds to the case where the element of Z N corresponding to the permutation index of c is 

0, which has order k = 1. For the cases where d ^ 0, continuing along the path and passing through q£ 
edges of the inverse image of c, we reach the node b ilP >, where p' = j + qd mod N. Clearly, we will be 
back to the starting node for the first time when qd = mod N, for the smallest value of q. (In that case 
by continuing the path we will just go over the same cycle of length q£.) By definition, q is the order of 
d in Z N , and thus q = k. Since the order of any element of a finite group divides the order of the group, 
s = N/k is an integer. It can then be easily seen that by starting from any of the nodes b^i, . . . , bi lS , we 
can partition the inverse image of c into s cycles of length k£ each. This completes the proof. ■ 

In what follows, we refer to the value k in Theorem Q] as the order of cycle c, and use the notation 
0(c) to denote it. 

Corollary 1: Consider the cyclic iV-lifting G of the Tanner graph G. Suppose that c is a cycle of length 
£ in G. The inverse image of c in G is the union of non-overlapping cycles, each strictly longer than £ 
iff 0(c) > 1; or equivalently, iff the permutation index of c, given in ©, is nonzero. 

B. Intentional Edge Swapping (IES) Algorithm 

Suppose that T is the set of all dominant trapping sets, and C(T) is the set of all the cycles in T. We 
also use the notations t and C(t) for a trapping set and its constituent cycles, respectively. For an edge e, 
we use C e (t) to denote the set of cycles in the trapping set t that include e. In the previous subsection, 
we proved that a cycle c in the base Tanner graph G is mapped to the union of larger cycles in the 
cyclically lifted graph G iff 0(c) > 1. To eliminate the dominant trapping sets, we are thus interested 
in assigning the edge permutation indices to the edges of C(T) such that 0(c) > 1 for every cycle 
c G C(T). To achieve this, we order the trapping sets in accordance with the increasing order of their 
critical number. We still denote this ordered set by T with a slight abuse of notation. Note that T may 
now include trapping sets with critical numbers larger than the minimum one. We then go through the 
trapping sets in T one at a time and identify and list all the cycles involved in each trapping set in C(T), 

1. e., C(T) = {c G C(t),Vt G T}. Note that C(T) is also partially ordered based on the corresponding 
ordering of the trapping sets in T. 

Example 1: Three typical trapping sets for Gallager A/B algorithms are shown in Fig. 1 Q. The (4,4) 
and (5, 3) trapping sets include one and three cycles of length 8, respectively, while the (4, 2) trapping 



8 



set has 2 cycles of length 6 and one cycle of length 8. 




Fig. 1. a) (5, 3) Trapping set b) (4, 2) Trapping set c) (4, 4) Trapping set. o = Variable Node, □ = Even degree Check Node and 
Odd degree Check Node. 



The next step is to choose proper edges of each trapping set to be swapped, i.e., to choose the edges to 
which nonzero permutation indices are assigned. In general, the policy is to select the minimum number 
of edges that can result in 0(c) > 1 for every cycle c in the trapping set t under consideration. 

Example 2: Going back to Fig. 1, for the (4, 4) trapping set, it would be enough to just pick one of the 
edges of the cycle of length 8 to eliminate this cycle in the lifted graph. For the (5, 3) and (4, 2) trapping 
sets, however, at least two edges should be selected for the elimination of all the cycles. A proper choice 
would be to select one edge from the diagonal and the other edge from one of the sides. 
Related to the edge selection, is the next step of permutation index assignment to the selected edges such 
that 0(c) > 1 for every cycle c G C(t). In general, we would like to have larger orders for the cycles. This 
in turn would result in larger cycles in the lifted graph. To limit the complexity, however, we approach 
this problem in a greedy fashion and with the main goal of just eliminating all the cycles in C(t), i.e., for 
each selected edge e, we choose the permutation index such that all the cycles C e (t) have orders larger 
than one. This can be performed by sequentially testing the values in the set Zi^^tO As soon as such an 
index is found, we assign it to e and move to the next selected edge and repeat the same process. 

2 More complex search algorithms with the goal of increasing the order of cycles may be devised. In this work however, no attempt has 
been made in this direction. 
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We call the proposed algorithm intentional edge swapping (IES) to distinguish it from "random edge 
swapping," commonly used to construct lifted codes and graphs. The pseudocode of the algorithm is 
given as Algorithm 1. At the output of Algorithm 1, we have the sets SwappedSet and IndexSet, which 
contain the edges of the Tanner graph that should be swapped, and their corresponding permutation indices, 
respectively. 

Algorithm 1 Intentional Edge Swapping (IES) Algorithm 

1) Initialization: Create the ordered sets T and C(T). Select N. ProcessedSet = 0, SwappedSet = 0, 
IndexSet = 0. 

2) Select the next trapping set t E T. 

3) Currents et = edges of C(t). 

4) CandidateSet = CurrentSet \ ProcessedSet. 

5) If CandidateSet = 0, go to Step 8. 

6) Select the edges £ from CandidateSet that should be swapped, and assign their permutation indices 
X from Zi^ N such that O(c) > 1 for every cycle c in C e (t). 

7) SwappedSet = SwappedSet U S, IndexSet = IndexSet U X, and ProcessedSet = ProcessedSet U 
CurrentSet. Go to Step 12. 

8) CandidateSet = CurrentSet \ SwappedSet. If CandidateSet = 0, Stop. 

9) Select an edge e from CandidateSet and assign a permutation index % £ Z\-±n to it such that for all 
cycles c £ C e (ProcessedSet U t), we have 0(c) > 1. If this is not feasible, go to Step 11. 

10) SwappedSet = SwappedSetUe, IndexSet = IndexSetUi, and CurrentSet = CurrentSet\C e (t). 
If CurrentSet ^ 0, go to Step 8. Otherwise, go to Step 12. 

11) CandidateSet = CandidateSet \ e, If CandidateSet ^ 0, go to Step 9. Else, stop. 

12) If all the trapping sets in T are processed, stop. Otherwise, go to Step 2. 

In Algorithm [H the search for edges to be swapped and the permutation index assignment to these 
edges are performed in two phases. The first phase is in Steps 4-6, where any edge from previously 
processed trapping sets is removed from the set of candidates for swapping. If the first phase fails, in 
that no edge exists as a candidate for swapping (CandidateSet = 0), then the algorithm switches to the 
second phase in Steps 8 - 9, where only previously swapped edges are removed from the candidate set 
for swapping. 

The process of permutation index assignment in Algorithm Q] involves the satisfaction of inequalities 
d 7^ for certain cycles, where d is given in (Q. In general, this is easier to achieve if the variables 
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involved in © are selected from a larger alphabet space. In fact, by increasing N, one can eliminate 
more trapping sets and achieve a better performance in the error floor region. 



C. Minimum Distance and Rate of Cyclic Liftings 

Consider an LDPC code C with anmxn parity-check matrix H. To prove our results on the minimum 
distance and the rate of a cyclic iV-lifting C of C, we consider an alternate parity-check matrix H' of C 
obtained by permutations of rows and columns of matrix H introduced in Subsection III-BL as follows: 



H' 



Aq An-i An-2 
A\ Aq An-\ 



A 2 A% 
A 3 A 2 



\ 



Ao An-i 
At A J 



(4) 



{i,j)th entry of Ad 



(5) 



An-2 An-3 An-4 

\ An-i An-2 Ans 
In ©, all the sub matrices Ad, d—0,...,N — l, have size m x n, and are given by 

1 if d = dij 
otherwise, 

where d^ is the permutation index corresponding to hij. The parity-check matrix H' is block circulant 
with the property that 

A0 + A1 + ... + An-i = H. (6) 

Theorem 2: If code C has rate r, then the rate r^ N ' of a cyclic iV-lifting C of C satisfies r^ N > < r, for 
N = 2*, g = l,2,3,.... 

Proof: Due to the block circulant structure of H', it can be written as 

M Af 
v J\f M 

where matrices M. and N are given by 

^ An—i ■ ■ ■ An_ +2 An_ +1 ^ 



H' 



(7) 



M 



Ai A ... ^^ +3 Ai^ +2 



2 



An 

2 



A x 



/ N N 
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and 



Af 



An 

2 



At 



Ai 



An 

2 



A 2 
A 3 



A 1 
A 2 



\ 



{An-, A N - 2 ... An +3 An +2 J fm>< ^ 
respectively. (Note that all indices 2 of ^ should be interpreted as modulo N.) Adding the second block 
column of © to the first, followed by adding the first block row to the second, we have 
M Af \ I M+Af AT \ f M+Af Af 
Af M J \Af + M M J \ M+Af 
For N = 2, A4 +Af = H, and since the rank of the matrix in ©, and thus the rank of H', is at least twice 
the rank of H, we have < r, and the proof is complete. For N > 2, it is easy to see that A4 + Af 
is also block circulant and can in turn be partitioned into four sub matrices, each of size jin x —n, as 



(8) 



follows: 

' M' Af' 
Af' M 

Replacing this in the rightmost matrix of ®, and performing similar block operations as in ([8]), we obtain 



M+Af 



( 



M! 


Af' 


Af 






Af' 


M' 








M! Af' 













Af' M! 



\ 



-+ 



\ 



M! + Af' Af' 


\ 




V 


M' + Af' 






M! + Af' Af' 









M' + Af' j 



(9) 



For N = 4, M' + Af' = H, and as the rank of the matrix in © is at least four times the rank of H, we 
have r( 4 ) < r, and the proof is complete. For N = 2 q > A, the same process of block column and row 
operations is repeated q times resulting in a block upper triangular matrix with the following structure 

/ H S 1)2 ... B ljJV _! B 1>N \ 
H ... B2N-1 B2N 





V 








H 




&N-1,N 

H 



(10) 



NmxNn 
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As the rank of the above matrix, and thus the rank of H', is at least N times the rank of H, we have 

r (N) < r j 

Corollary 2: If matrix iJ has full rank, then r( N > = r, for iV = 2 g , q = 1, 2, 3, . . .. 
Proof: If has full rank, then the matrix in (flOl) is also full-rank, and so is H'. This implies = r. 

■ 

It is easy to find counter examples to demonstrate that Theorem [2] and Corollary [2] do not always hold 
for odd values of N or even values that are not integer powers of two. 

Theorem 3: If code C has minimum distance d min , then the minimum distance d^ n of a cyclic iV-lifting 
C of C satisfies d min <d [ ^ n <N d min , for N = 2 q , g = 1, 2, 3, . . .. 

Proof: We first prove the lower bound. Consider a codeword with minimum Hamming weight 
^min = Wft(c^) in C. Let c^) = (c , c 1; . . . , c Ar _ 1 ), where the subvectors c i; i = 0, 1, . . . , iV — 1, are of 
size n each. Based on H'c/ N ^ = 0, we have 

N-l 

EAv-^mod^ = Q, J = 0,...,iV-l. (11) 

Adding all the equations in ([111) for different values of j, and exchanging the order of summations over 
i and j, we obtain 

N-l (N-l \ JV-1 /JV-l \ 

E E ^-i+imodiv & = E E^U=Q- (12) 

i=0 \J=0 / i=0 \k=0 J 

Based on (J6]), this implies that c = Cq+Ci + - ■ •+c 7V _ 1 is a codeword of C. Moreover, Wh(c) < Wft(c^) = 
d^. If c 7^ , this means d min < d^ n . 

If c = and N = 2, then c = This along with H'cW = results in Hc = 0, and thus c G C. 
Since c 7^ 0, this implies <i m j n < >V/ l (c ) = d^ in /2 < d^ in and the proof is complete. 

If c = and N > 2, through a number of steps, we demonstrate that either there exists a subset s of 
Zq-^n-i, such that the vector y = 2~2 ies c,i is nonzero and is in C, or = c i+Ar / 2 , fori = 0, 1, . . . , N/2 — 1. 
For the former case, d min < Wh(y_) < V\4(c s ) < W/^c^) = d^ n , where c s is defined as the vector 
of size \s\ • n obtained by the concatenation of vectors c i , i E s. This proves the lower bound. For the 
latter case, the problem is reduced to that of a lifting with degree N/2. Iterating the same process, we 
either find a nonzero vector of C with Hamming weight less than d^ n or reach to a point where all the 
constituent vectors c ,c 1; . . . ,c N _ x of c are equal. In this case, vector y = c = c 1 = ■ • ■ = c N _ x is 
nonzero and is in C. We thus have d min < W h {y_) = d^ n /N < d^ n . 

Here, we explain the first step for demonstrating that when c = and N > 2, either there exists a 
subset s of Z _>jv-i> such that the vector y = J2ies^i ls nonzero and is in C, or c t = c i+N , 2 , for % = 
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0, 1, ... , N/2 — 1 . The other steps are similar and omitted to avoid redundancy. (Note that the first step 
suffices to prove the claim for N = A. For larger values of N further steps are required.) Consider the 
vectors c-, i = 0, 1, . . . , N — 1, defined by 

i-l+N/2 



E %modjv ( 13 ) 



Using equations (fTTI) . it is easy to see that vector c'^ = (c^c^, . . . , Cj V _ 1 ) satisfies H'g?( N > = 0, and is 
therefore in C. Moreover, from c = , we have 

Pi = cj +JV/2 , for z = 0,...,iV/2-l. (14) 

Applying this to equation H'd^ N ' = 0, we obtain 

N/2-1 

(^N-i+j mod n + ^N/2-i+j mod n) -i = - ' i = 0' • • • ' — 1 • (15) 

i=0 

Adding the equations in (fl~5T) for different values of j and switching the summations with respect to i and 
j, we obtain 

'N/2-1 \ /JV-1 

E ^ E^) = -' (16) 

i=o y Vfc=o / 

which implies that the vector c' = c + c' x + h 9/.N/2-1 * s a codeword of C. On the other hand, using 

definition (PT3T) . we have c' = c + c 2 + • • • + c N _ 2 . (The subset s in this case is {0,2, . . . , N — 2}.) 
Thus W h {d) < Wh^Cz, . . .,c N _ 2 )) < W ft (cW), implying d min < if ^ Q. If c' = 0, then 
J2ies^-i = Q ° ver both the even and odd subsets s of Z _>jv_i. For N — A, this means c = c 2 and 
£i = £3- Replacing these in (fTTT) . we have 

U ^ + Al \ I Cn \ 

0. 

Now the problem is reduced to that of iV = 2, which means y = c + c 1 E C and either y ^ 0, or c = c v 
In both cases, the lower bound is proved. 

To prove the upper bound, consider a codeword c G C with W&(c) = d min . Define a vector of size 
n • iV by = (c, . . . , c). It is easy to see based on © that c^ satisfies H'cS N > = and is thus in the 
cyclic ^-lifting C of C. We therefore have d^ n < W h (c^) = N-W h (c) = N- d min . U 

It is important to note that Theorems 1 and 2 of ifTOll are special cases of Theorems [2] and [3] of this 
paper, respectively, where N = 2 q = 2. 



A 


+ A 2 


A! 


+ A 3 \ 




A, 


+ A 3 


A 


+ A 2 I 
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IV. NUMERICAL RESULTS 

In this section, we apply the IES algorithm of Subsection IIII-BI to three LDPC codes to eliminate 
their dominant trapping sets over the BSC. The codes are: the (155,64) Tanner code [181, a (504,252) 
randomly constructed regular code [fT2l . and an optimized (200, 100) randomly constructed irregular code. 

Example 3: For the (155, 64) Tanner code under Gallager B algorithm, the most dominant trapping set 
is the (5, 3) trapping set, shown in Fig. \T\ with critical number 3. We apply the IES algorithm to this code 
to design cyclic iV-liftings for N = 2,3,4, and 5. The FER curves of the designed codes are presented 
in Fig. [2] along with the FER of the base code. 

A careful inspection of Fig [2] shows that using a 2-lifting, the slope of the curve changes from 3 to 4, 
an indication that all (5, 3) trapping sets are eliminated. In this case, (4, 4) trapping sets play the dominant 
role. Further increase of N to 3 and then 4, only causes a downward shift of the curve (with no change 
of slope), an indication that the minimal critical number remains at 4 for the 2 lifted codes and increasing 
the degree of lifting just reduces the number of (4, 4) trapping sets. Increasing N to 5 however, eliminates 
all the (4,4) trapping sets and the slope of the FER curve further increases to 5. The dominant trapping 
sets for the 5-lifting are (5, 5) trapping sets. 

It is important to note that for N = 2, the performance of the designed code is practically identical to 
that of the code designed in Example 3 of ifTOl based on a 2-cover of the Tanner code. There are however 
no results reported in IfTOl for covers of larger degree. 

For comparison, we have also included in Fig |2l the FER of a random 5-lifting of the Tanner code. 
As can be seen, the error floor performance of this code is significantly worse than that of the designed 
5-lifting. In particular, the slope of the random lifting is just 4 versus 5 for the designed lifting. 

The code rates of the designed iV-liftings are: 0.4065, 0.4043, 0.4032, and 0.4026, for N = 2 to 5, 
respectively. The small decrease in the code rate by increasing the degree of lifting is a consequence of 
the fact that the original parity-check matrix of the Tanner code is not full rank. The rate of the Tanner 
code itself is 0.4129. 

It is also worth noting that while the girth of the iV-liftings, N = 2,3,4, remains the same as that of 
the Tanner code, i.e., g = 8, for the 5-lifting, the girth is increased to 10. 

Example 4: In this example, we consider a regular (504,252) code from [12] decoded by Gallager B 
algorithm. The dominant trapping sets in this case have critical number 3 and include (3,3), (4,2), and 
(5, 3) trapping sets among others. The IES algorithm is used to design cyclic iV-liftings of this code 
for N = 2 to 6. The FER results of the liftings and the base code are reported in Fig. [3] Again, the 
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Fig. 2. Comparison of the FER performance of the Tanner code and its liftings over the BSC (Example O. 

performance of the 2-lifting is similar to that of the code designed in [10]. All (3,3) trapping sets are 
eliminated in the 2-lifting, but the survival of other trapping sets with critical number 3 keeps the minimal 
critical number at 3, and thus no change of FER slope compared to the base code is attained. Increasing 
N to 3, however, eliminates all the trapping sets with critical number 3 and changes the slope of the FER 
to 4. The dominant trapping sets in this case are (4, 4) sets. Further increase of N to 4 and 5 only reduces 
the number of (4, 4) trapping sets and thus results in a downward shift of the FER curve. For N = 6, the 
algorithm can eliminate all the (4,4) trapping sets, and thus increases the slope of the FER curve to 5. 
The dominant trapping sets in this case are (5, 5) sets. 

For comparison, in Fig. |3l we have also shown the performance of a random 6-lifting of the (504, 252) 
code. As can be seen the performance of this code in the error floor region is far poorer than that of 
the designed 6-lifting. In particular, the slope of the FER curve for this code is only 3 versus 5 for the 
designed code. 

In this example, the parity-check matrix of the base code is full-rank, and all the liftings have the same 
rate of 0.5 as the base code. 

Noteworthy is that while the 2-lifting has the same girth of g = 6 as the base code, the girth for 
Af-liftings, N = 3 to 6, is increased respectively to 8, 8, 8 and 10. 

Example 5: In this example, we consider a randomly constructed rate- 1/2 irregular (200, 100) code 
as the base code. The degree distributions for this code, optimized for Gallager A algorithm over the 



- Tanner (155.64) 

- IES 2-lifting (310,126) 

- iES 3— lifting (465,188) 

- IES 4-lilting (620,250) 

- IES 5-lilting (755,312) 

- Random 5— lifting (755,312) 



Slope 3 
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Fig. 3. Comparison of the FER performance of the regular (504, 252) code and its liftings over the BSC (Example [4j. 

BSC [1J, are X(x) = 0.1115x 2 + 0.8885x 3 and p(x) = 0.26x 6 + 0.74x 7 . The code has g = 6. This code 
has a wide variety of dominant trapping sets under Gallager A algorithm, all with critical number 3. 

We apply the IES algorithm to this code to design a cyclic 13-lifting of length 2600, rate 0.5 and g = 6. 
The FER curves of the lifted code and the base code are presented in Fig. 0] As can be seen, the lifted 
code has a much better error floor performance compared to the base code. In fact, the minimum critical 
number for the lifted code is 5 versus 3 for the base code. For comparison, a rate- 1/2 code of block length 
2600 with the same degree distribution and g = 6 is constructed. The performance of this code is also 
given in Fig. HI Clearly the performance of the lifted code is significantly better in the error floor region. 
In particular, the slope of the FER curve for the random code is the same as the base code and much less 
than that of the lifted code. 

Example 6: It is known that codes designed for a certain channel/decoding algorithm would also 
perform well for other channel/decoding algorithms [7]. In this example, we show that cyclically lifted 
codes designed for Gallager B algorithm in Examples [3] and 0] also perform very well on the binary-input 
AWGN channel under min-sum algorithm. The FER results for the 5-lifting of the Tanner code and the 
6-lifting of the Mac Kay code are reported in Figures [5] and [6l respectively. In each figure, the performance 
of the corresponding base code and a similar random code (same block length and degree distributions) 
is also presented. One can see that at high SNR values, the designed codes perform far superior to the 
corresponding base codes and random codes. In particular, they show no sign of error floor for FER 
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Fig. 4. FER curves of the irregular (200, 100) code, its IES cyclic 13-lifting and a random irregular (2600, 1300) code with the same 
degree distribution over the BSC (Example |5). 

values down to about 1CT 8 , and their FER decreases at a much faster rate compared to the base codes 
and random codes. 



In this work, cyclic liftings are proposed to improve the error floor performance of LDPC codes. The 
liftings are designed to eliminate the dominant trapping sets of the code by eliminating their constituent 
short cycles. The design approach is universal in that it can be applied to any decoding algorithm over any 
channel, as long as the dominant trapping sets are known and available. In addition, the liftings have the 
same degree distribution as the base code and are implementation friendly due to their cyclic structure. 
For base codes with full-rank parity-check matrices, the liftings also have the same rate as the base code 
and the performance improvement is achieved at the expense of larger block length. Compared to random 
codes or random liftings with the same block length and degree distribution, the designed codes perform 
significantly better in the error floor region. 

While the cyclic liftings in this work were designed for Gallager A/B algorithms over the BSC, they 
also performed very well over the BIAWGN channel. In particular, the designed codes substantially 
outperformed similar random codes in the high SNR region. 



V. CONCLUSIONS 
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Fig. 5. FER performance of the Tanner code, its IES cyclic 5-lifting and a random (775, 310) code with the same degree distribution over 
the BIAWGN channel (Example [6} . 




Eb/No (dB) 



Fig. 6. FER performance of the regular (504, 252) code, its IES cyclic 6-lifting and a random (3024, 1512) code with the same degree 
distribution over the BIAWGN channel (Example^. 
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